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Abstract 

We report measurements of time dependent decay rates for B'^{B'^) Z)*-*)^7r^ decays and 
extraction of CP violation parameters containing ^3. Using fully reconstructed D^*^Tr events from 
a 140 fb~^ data sample collected at the T(45) resonance, we obtain the CP violation parameters for 
D*Tr and Dtt decays, 2i?^(«)^ sin(2^i + ^3^: Sd(*)^), where Rjji*)^ is the ratio of the magnitudes of 
the doubly-Cabibbo-suppressed and Cabibbo-favoured amplitudes, and Sj^(t)^ is the strong phase 
difference between them. Under the assumption of Sjji*)^^ being close to either or 180°, we 
obtain |2i2D,^ sin(2</>i + ^3)! = 0.060 ± 0.040(stat) ± 0.019(sys) and |2ii£,7r sin(2^i + ^3)! = 0.061 ± 
0.037(stat) ± 0.018(sys). 

PACS numbers: 12.15.Hh,13.25.Hw 
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The good agreement between direct measurements of sin 2(/)i P, 0] and the outcome of 
global fits to the CKM quark mixing matrix elements Q strongly supports the standard 
model explanation of CP violation. To determine whether it is the complete description or 
whether additional factors come into play, further measurements of other CKM parameters 
are required. Among these parameters 03 is of particular importance. The measurements 
of time-dependent decay rates of B^{B^) — > D(*^^7r^ provide a theoretically clean method 
for extracting sin(20i + 03), since loop diagrams do not contribute to these decays 

There are two ways for a state which is initially 5° to be found as D^*'^~'k^ at a later time 
t. It can occur either directly through a Cabbibo-favoured decay (CFD) or through mixing 
followed by doubly- Cabbibo-suppressed decay (DCSD), as shown in Fig. Q Interference of 
the two processes introduces the term containing 03 to the time dependent decay rates, 
which are given by |5i, ISI] 



. u 



It , D' 



u 



d , d - B ^ ^ It 

(a) CFD (b) DCSD 

FIG. 1: Contributions to 5° D^*^ 'TT~^ can come either (a) from CFD or (b) from mixing 
followed by DCSD. 



P(E° L)(*)+7r-) = c [1 - cos Amt - 23p sin Amt] 
P{B^ ^ D^*^-Tx+) = c[l + cosAmt + 2^psmAmt] 
P(i?° ^ dW+tt") = c[l + cosAmt + 2$5psinAmt] 
P(fi°^DW-7r+) = c[l - cosAmt -2'^psmAmt] 

(1) 

where c = (6^*/"^^° )/2tbo with r^o denoting the lifetime of the neutral B meson and Am is the 
B^-B^ mixing parameter. The p and p are defined as p = {q/p){A{B^ D^*^^n^) / A{B^ — >■ 
P)(*)-7r+)) and p = (p/g)(^(P° ^ D^*^+7i-)/A{B^ DW+tt")), where p and q relate the 
mass eigenstates to the flavour eigenstates in the neutral B meson system jsj. They lead to 
CP violating terms Q'p = — (— l)^i?sin(20i+03— 5) and 55p = (—l)'^i?sin(20i-|-03-|-(5), where 
R is the ratio of the magnitudes of the DCSD and CFD amplitudes (here the magnitudes 
of both the CFD and DCSD amplitudes are assumed to be same for B^ and B^ decays), S 
is the strong phase difference between DCSD and CFD, and L is the angular momentum of 
the final state (1 for D*n and for Dtt). R and 6 are not necessarily the same for D*n and 
Dn final states, and are denoted with subscripts, D*7c and Dtc, in what follows. 

This study uses a 140 fb~^ data sample, which contains 152 million BB events, collected 
with the Belle detector p] at the KEKB collider The selection of hadronic events is 
described elsewhere 0. 

For the B^ D*~^n~ event selection, we use the decay chain D*^ D^n^, and — >■ 
K~TT'^, K^Ti^TT^ or K~7i^7i^'ri^ (charge conjugate modes are implied throughout this paper). 
For the B^ D^tt~ event selection, we use — > K~7c^tt^ decays. Charged tracks except 
the slow vr^ in the D*^ — > D^vr"'" decay are required to have a minimum of one hit (two hits) 
in the r-0 (z) plane of the vertex detector in order to allow precise vertex determination. To 
separate kaons from pions, we form a likelihood for each track, CK{n)- The kaon likelihood 
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ratio, P{K/n) = Ck/{Ck + -Cvr), has values between (likely to be a pion) and 1 (likely 
to be a kaon). We require charged kaons to satisfy P{K/7r) > 0.3. No such requirement is 
imposed to select charged pions coming from D decays. 

For selection, the invariant mass of the daughter particles is required to be within 
±16.5MeV/c2, ±24.0MeV/c2, and ilS.SMeV/c^ of the nominal mass, for K-7r+, 
K^tt^tt'^, and K^TT~^7T~^'n'~ modes, respectively. These intervals correspond to ±3(T, where a 
is the Monte Carlo determined invariant mass resolution. For the D~^, the invariant mass is 
required to be within ±12.5 MeV/c^ of the nominal mass. For the D° K' '7i~^n^ recon- 
struction, we further require the 7r° momentum to be greater than 200 MeV/c in the T(45') 
rest frame, and the ratio of the second to zeroth Fox- Wolfram moments [10| R2 to be less 
than 0.55. We require R2 < 0.5 for —>■ K'tt^tt^ . We use a mass- and vertex-constrained 
fit for and a vertex-constrained fit for . 

The D*^ is reconstructed by combining candidates with a slow 7r+. Here, slow pions 
are required to have momentum less than 300MeV/c in the T(4S') rest frame. The D* 
candidates are required to have a mass difference AM = Mdo^^ — M^o within ±7MeV/c^, 
±2MeV/c^, or ±4MeV/c^ of the nominal value, for the K~7r~^, K~it~^tc^, and /^"Tr+Tr+Tr" 
modes respectively. 

We reconstruct B candidates by combining the candidate with a tt^ candidate 

satisfying P{K/tx) < 0.8. We identify B decays based on requirements on the energy differ- 
ence /S.E ^J2iEi — -Ebeam and the beam-energy constrained mass Mbc = i/-Ebeam ~ (SiPi)^; 
where -Ebeam is the beam energy, pi and Ei are the momenta and energies of the daugh- 
ters of the reconstructed B meson candidate, all in the T(45') rest frame. If more than 
one B candidate is found in the same event, we select the one with best D vertex quality. 
We define a signal region in the AE-M^c plane of 5.27GeV/c^ < Mbc < 5.29GeV/c^ and 
\A.E\ < 0.045 GeV, corresponding to about ±3o" of both quantities. For the determination 
of background parameters, we use events in a sideband region defined by Mbc > 5.2 GeV/c^ 
and —0.14 GeV < AE < 0.20 GeV, excluding the signal region. 

Charged leptons, pions, and kaons that are not associated with the reconstructed D^*^ti 
decays are used to identify the fiavour of the accompanying B meson. The algorithm ^ 
leads to two parameters, q and r, where g = +1 indicates b hence B^ and g = — 1 indicates 
b hence B^. The parameter r is an event-by-event dilution factor ranging from r = for no 
fiavour discrimination to r = 1 for unambiguous fiavour assignment. More than 99.5% of 
the events are assigned non-zero values of r. 

The decay vertices of the B — >■ D^*'^7i are fitted using the momentum vectors of the D and 
TT (except the slow vr from D* decay) and a requirement that they are consistent with the 
interaction region profile. For the decay vertices of the tagging B meson, the remaining well 
reconstructed tracks in the event are used. Tracks that are consistent with Kg decay are 
rejected. The proper-time difference between the fully reconstructed and the associated B 
decay is calculated as At = {zj-ec — 2;tag)/c/37, where z^ec and Zt^g are the z coordinates of the 
two B decay vertices and jSj = 0.425 is the Lorentz boost factor at KEKB. After application 
of the event selection criteria and the requirement that both i?'s have well defined vertices 
and |At| < 70 ps (~ A5tbo), 7763 and 9351 events remain as the D*7t and Dn candidates, 
respectively. The signal fractions of the samples, which vary for different r bins, are 96% 
for D*7i and 91% for Dti. 

Unbinned maximum likelihood fits to the four time dependent decay rates are performed 
to extract and Qp. We minimize — 2X)ilnLi where the likelihood for the i-th event is 
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given by 



Li — (l — fol) [/sig-Psig ® -Rsig + (1 — /sig)-Pbkg ® -Rbkg] 
+ /ol-Pol- 

The signal fraction /gig is determined from the {AE, Mbc) value of each event. The signal 
distribution is the product of the sum of two Gaussian in AE and a Gaussian in Mbc; 
that for the background is the product of a first order polynomial in AE and an ARGUS 
function in Mbc- 

The At distribution is modeled by a core distribution convolved with resolutions. A small 
number of events have poorly reconstructed vertices resulting in a very broad distribution At. 
We account for the contributions from these "outliers" by adding a Gaussian component PqI 
with a width and fraction determined from the B lifetime analysis 3|. The At resolution, 



denoted by i?sig and i?bkg for the signal and background, is determined on an event-by-event 
basis, using the estimated uncertainties on the z vertex positions [l3| . 
The signal At distributions are given by 

Psig(g = -l,Z}W±7r^) 

= (1 - w_)P{B° D(*)±7r^) + w+P{BO DW^tt^) 
= (1 - w+)P{BO D(*)±7r^) + w^P{B'^ DW^tt^) 

(2) 

where W7_ and w+ are wrong tag fractions for the q = —1 and q = +1 samples, respectively. 
P's are given by Eq. Qwith t and c replaced by At and (e~''^*'/'^so)/4r^o, respectively. 

The background At distribution is parameterized as a sum of a 5-function component 
and an exponential component with a lifetime Tbkg 

^'bkg = fL,m - pi,^) + iL_^e-|^*---J/^-^ (3) 

bkg 

where /bkg is the fraction of events contained in the 5-function, and /ibkg and /ibkg are the 
mean values of |At| in the ^-function and exponential components, respectively. 

While the tagging side should have no asymmetry if the flavour is tagged by primary 
leptons, it is possible to introduce a small asymmetry when daughter particles of hadronic 
decays such as D^*^7i are used for the flavour tagging, due to the same CP violating effect, 
which is the subject of this paper 0]. This effect is taken into account by replacing the 
coefficients of sin Amt in Eqs. [Uby Qp — '^p', Qp — Qp', Qp — Qp', and Qp — Qp', respectively. 
Here the and Qp' represent the CP violating effect due to the presence of B^ DX 
and B^ — »• DX amplitudes in the flavour tagging side. Note that unlike the Qp and Qp, 
which are rigorously defined in terms of B^ D^*^^^^ and B^ — >■ D*^*)^7r'^ amplitudes, ^p' 
and are effective quantities that include effects of the fraction of P ^ DX components 
in the tagging B decays and all experimental effects of subsequent behaviour of D mesons. 
Therefore, these quantities must be determined experimentally. 

The values of '^p' and S^p' are determined in each of six r bins by fitting the At dis- 
tributions of a D*lu control sample using the signal distributions of Eq. El and set- 
ting 53p and 53p to zero. Since the D*lu final states have specific fiavour, any observable 



6 



asymmetry must originate from the tagging side. The results for the combined r bins are 
2 V = 0.038 ± 0.014(stat) ± 0.005(sys) and 2Qp' = 0.002 ± 0.014(stat) ± 0.009(sys). 

The procedures for At determination and flavour tagging are tested by extracting r^o and 
Am. When all four signal categories in Eq.^are combined, the signal At distribution reduces 
to an exponential lifetime distribution. We obtain r^o = 1.583 ± 0.029 ps (1.575 ± 0.032 ps) 
for the D*7i (Dtt) samples, in good agreement with the world average (1.542 ±0.016 ps) j3[. 
Combining the two CFD-dominant modes and the two mixing-dominant modes and ignoring 
the CP violating terms, the asymmetry behaves as cos AmAt. We obtain Am = 0.490 ± 
0.015 ps~^ (0.483 ± 0.014 ps~^) for the 0*71 (Dn) samples, also in good agreement with the 
world average (0.489 ± 0.008 ps~^) 3]. The same fits also provide wrong tag fractions W- 
and w+ in each r bin for both D*7C and Dtt data samples. The errors of our results are 
statistical only. 

We then perform fits to determine the and by fixing r^o and Am^ to the world 
average values and using w_, Qp', and Qp' for each r bin, as obtained from the above fits. 
The resuhs are 2QpD*n = 0.011 ±0.057, 2'^pd*^ = -0.109±0.057, 2'^pDn = -0.037±0.052, 
and 2S5p£)7r = 0.087 ± 0.054. The errors are statistical only. The At distributions for the 
subsamples having the best quality flavour tagging (0.875 < r < 1.000) are shown in Fig. |21 
for the D*TT and in Fig. El for the Dtc samples, respectively. 

The systematic errors come from i) the uncertainties of parameters which are constrained 
in the fit, including At resolution parameters, background parameters, wrong tag fractions, 
and physics parameters; ii) uncertainties of the tagging side asymmetries; iii) fit biases 
induced by the vertexing and other unknown factors. For item i) , we repeat the fits varying 
each parameter value by ±la. To estimate item ii), we repeat the fits by varying the Qp' and 
Qp' by their errors. Errors are not explicitly assigned for item iii), since they are included 
in the errors of '^p' and Qp' from the D*/z/ control sample fit (item ii). Table |l] summarizes 
the systematic errors. 



TABLE I: Systematic errors in the 2R s'm{2(j)i ± (/)3 ± 5) extractions. 



Sources D*^ Dtt 

Signal At resolution 0.014 0.013 

Background At shape 0.001 0.003 

Background fraction 0.002 0.001 

Wrong tag fraction 0.006 0.006 

Vertexing 0.005 0.005 
Physics parameters (Am,r^o) 0.001 0.002 

Tagging side asymmetry 0.009 0.009 

Combined 0.019 0.018 



2RD*n sin(20i ± (^3 ± Sd*^) 
2Rd*^ sin(20i ± ^3 - Sd*^) 
2Rn^ sin(20i ± 03 ± ^dtt) 
2i?D7r sin(20i ± 03 - ^Dtt) 



= 0.109 ±0.057 ±0.019, 
= 0.011 ±0.057 ±0.019, 
= 0.087 ±0.054 ±0.018, 
= 0.037 ±0.052 ±0.018. 
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FIG. 2: At distributions for the D*it data in the 0.875 < r < 1.000 flavour tagging quahty bin. 
(a) B° ^ D*+n-, (b) D*-Tr+, (c) B° D*+Tr-, (d) Z)*-7r+. Curves show the fit 

results with the entire event sample, hatched regions indicate the backgrounds. 




FIG. 3: At distributions for the Dtt events in the 0.875 < r < 1.000 flavour tagging quality bin. 
(a) B° D+TT-, (b) B° D-TT+, (c) B° D+tt-, (d) B° D--K+. Curves show the fit results 
with the entire event sample, hatched regions indicate the backgrounds. 



The first and second errors are statistical and systematic. At present, the statistical errors 
are too large to allow any meaningful conclusion to be drawn. However, it is interesting 
to consider how the four results can be combined using knowledge of R and 5 to irnprove 
the precision of sin(20i + ^3). Several methods have been proposed to measure i? A 
method that compares the branching fractions of 5° DI*^^-k^ and B'^ D^*^~tt~^ and 
uses factorization relation gives Rd^ = 0.0237 ± 0.0050 and Ro^n = 0.0180 ± 0.00670. 
The present errors are too large to conclude that the two R values are equal. On the other 
hand, there are solid theoretical grounds for assuming 60*^ and 6d-„ to be very small and 
therefore equal jl3]. However, some argue that there is an ambiguity of 180° between 6d*tt 
and ^Dn fs*]. Assuming is close to either 0° or 180°, we obtain |2i?£).7r sin(20i + 03)| = 

0.060 ±0.040(stat)±0.019(sys) and |2i?D^ sin(20i + ^s)] = 0.061 ± 0.037(stat) ± 0.018(sys). 



8 



In summary, we measure the time dependent CP violation parameters 2i?sin(20i + 03±5) 
for the 5° (5°) ^ D^*^^7i^ decays using 152 milhon BB events. Under the assumption of 
6oM^ being close to either 0° or 180°, we obtain |2i?D*^ sin(20i + 03)| = 0.060±0.040(stat)± 
0.019(sys) and |2i?^^ sin(20i + 03)| = 0.061 ± 0.037(stat) ± 0.018(sys). 
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